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Abstract 

Background: Tomato {Solanum lycopersicum) establishes a beneficial symbiosis with arbuscular mycorrhizal (AM) 
fungi. The formation of the mycorrhizal association in the roots leads to plant-wide modulation of gene expression. 
To understand the systemic effect of the fungal symbiosis on the tomato fruit, we used RNA-Seq to perform global 
transcriptome profiling on Moneymaker tomato fruits at the turning ripening stage. 

Results: Fruits were collected at 55 days after flowering, from plants colonized with Funneliformis mosseae and from 
control plants, which were fertilized to avoid responses related to nutrient deficiency. Transcriptome analysis 
identified 712 genes that are differentially expressed in fruits from mycorrhizal and control plants. Gene Ontology 
(GO) enrichment analysis of these genes showed 81 overrepresented functional GO classes. Up-regulated GO classes 
include photosynthesis, stress response, transport, amino acid synthesis and carbohydrate metabolism functions, 
suggesting a general impact of fungal symbiosis on primary metabolisms and, particularly, on mineral nutrition. 
Down-regulated GO classes include cell wall, metabolism and ethylene response pathways. Ouantitative RT-PCR 
validated the RNA-Seq results for 12 genes out of 14 when tested at three fruit ripening stages, mature green, 
breaker and turning. Quantification of fruit nutraceutical and mineral contents produced values consistent with the 
expression changes observed by RNA-Seq analysis. 

Conclusions: This RNA-Seq profiling produced a novel data set that explores the intersection of mycorrhization 
and fruit development. We found that the fruits of mycorrhizal plants show two transcriptomic "signatures": genes 
characteristic of a climacteric fleshy fruit, and genes characteristic of mycorrhizal status, like phosphate and sulphate 
transporters. Moreover, mycorrhizal plants under low nutrient conditions produce fruits with a nutrient content 
similar to those from non-mycorrhizal plants under high nutrient conditions, indicating that AM fungi can help 
replace exogenous fertilizer for fruit crops. 
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Background reciprocal reward for the best partner [3]. Furthermore, 

Tomato {Solanum lycopersicum), like other solanaceous research on tomato mycorrhizal roots has identified 

species, establishes a beneficial root symbiosis with plant genes as functional markers for AM symbiosis and 

arbuscular mycorrhizal (AM) fungi of the Glomeromy- shown that they are expressed in different cell types 

cota [1]. The plant roots accommodate the AM fungi, [4,5]. In addition to the expected effects on roots, 

which improve nutrient uptake of their plant hosts, and, mycorrhization also affects gene expression in tomato 

in turn, receive fixed carbon from the plant [2], with a shoots and leaves, modulating genes involved in diverse 

metabolic processes such as defence, transport and hor- 
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Tomato represents an excellent model plant for under- 
standing the long-distance effects of AM fungi on organs 
other than roots, as it exhibits substantial increases in 
growth and fruit yield in response to mycorrhization 
[9-11]. In addition, extracts of tomato fruit produced 
from mycorrhizal plants showed an anti-oestrogenic ac- 
tivity likely attributable to their higher contents of lyco- 
pene [12]. 

Analysis of the tomato {Solanum lycopersicum) gen- 
ome [13] has showed that many tomato genes are specif- 
ically related to the evolution of fleshy fruits, and that 
intensive breeding has selected some agriculturally rele- 
vant traits. Recent work has characterized the molecular 
mechanisms that govern the transition of tomato fruit 
from a green to a fully red ripe fruit, via the intermedi- 
ate stages termed breaker, turning and pink. Tomato rip- 
ening involves increased levels of aromatic compounds 
and carotenoid pigments, softening of the cell wall and 
decreased photosynthetic activity. Ethylene, the key hor- 
mone regulating ripening in climacteric fruit, controls 
most of these changes through a complex mechanism 
[14,15]. Accumulating data have helped to define the 
roles played by specific genes [16] in specific fruit cell 
types [17], and also have extended the analysis to tomato 
epigenetic and proteomic dynamics [18]. 

To better explore the intersection between mycorrhiza- 
tion and fruit development in tomato, we used RNA-Seq 
to profile transcription changes [19] in fruits originating 
from mycorrhizal plants, a condition which has so far not 
been analysed in a crop plant of high agricultural interest. 
We found that the fruits of mycorrhizal plants show regu- 
lation of transcripts characteristic of climacteric fleshy 
fruits and a core set of genes typically modulated by sym- 
biosis. Moreover, our quantitative analyses showed that 
fruits from mycorrhizal and non-mycorrhizal fertilized 
plants have comparable levels of carotenoids, phenols, and 
minerals (i.e. phosphorus, sulfur and potassium) indicating 
that the AM symbiosis benefits fruit quality. Thus, AM 
fungi can serve as surrogates for fertilization, providing 
similar benefits at the root and fruit levels. 

Results 

Plant sampling for RNA-Seq 

To examine the transcriptome of tomato fruits from 
mycorrhizal plants, we collected fruits from plants colo- 
nized by the AM fungus Funneliformis mosseae (previ- 
ously named Glomus mosseae) and from non-inoculated 
plants. AM fungi enhance nutrient uptake by their hosts 
[20]; to control for nutrition-dependent effects, we pro- 
vided increased nutrients to the non-inoculated plants, 
to produce two groups of plants that have a similar 
phenology and nutritional status. Thus, we watered the 
mycorrhizal plants (MYC plants) with a low-phosphorus 
(P), nitrogen (N) and sulfur (S) solution, which allowed 



high mycorrhization levels and healthy plant growth 
[8,21] but we watered the control, non-mycorrhizal 
plants (CONT plants) with higher concentrations of P, S 
and N. The MYC and CONT tomato plants grown 
under these different conditions showed similar growth 
(Figure 1) and similar P, S and K contents in their fruits 
(Additional file 1: Table SI), suggesting that the my- 
corrhizal fungus successfully compensated for the low 
fertilization treatment. We also measured some pheno- 
logical features of the plant, including fruit yield, and 
found that the MYC and CONT fruits showed compar- 
able values (Additional file 1: Table S2). The one excep- 
tion was fruiting time, which was longer for mycorrhizal 
plants, leading to the production of more fruit. By con- 
trast, plants that were not inoculated with F. mosseae 
but were fertilized with a similar, low-P and -N regimen 
to MYC plants were strongly affected by the nutrient 
limitations, showing impaired growth (Figure 1) and de- 
layed fruit production, leading to the production of 
fewer fruit. 

To assess the effect of mycorrhization on a key stage 
of ripening, we collected fruits at the turning stage, 
when the green color starts to change on 10 to 30 per- 
cent of the fruit surface, indicating the beginning of ca- 
rotenoid biosynthesis [22]. Because the MYC and CONT 
plants have comparable vegetative periods, we harvested 
fruits from both conditions at 55 days after flowering. 
The MYC plants showed a colonization frequency (F) of 
65% and a mycorrhizal intensity (M) of 36% in the root 
system (Additional file 2: Figure SI). We observed no 
traces of fungal colonization in the CONT roots. Two 
biological replicates for each condition were collected, 
following the current RNA-Seq standards (http://enco- 
deproject.org/ENCODE/protocols/dataStandards/ENCO- 
DE_RNAseq_Standards_V1.0.pdf). The samples were 
then used for deep single-end sequencing (read length: 
100 bp). 

Analyses of RNA-Seq data: read number, transcriptome 
coverage and total expressed genes 

To identify differentially expressed transcripts, we first fil- 
tered the sequences, estimated the transcriptome cover- 
age, and identified expressed genes. The raw reads (100 
bases, single-end) obtained from lUumina HIseq (FAS- 
TERIS Co. Geneva, Switzerland) were filtered (lUumina 
passed-filter call) and further checked for sequence con- 
taminants with the fastQC application. Contaminant-free, 
filtered reads ranging from 11 to 16 million for each sam- 
ple (Table 1) were mapped with Bowtie/TopHat to the to- 
mato genome sequence SL2.40 (ITAG2.3) [13]. Raw read 
counts were obtained from alignment files by counting 
with HTSeq software. Based on the sum of transcript 
lengths, as reported in the current SL2.40 annotation 
(41,393,518 bp; which does not currently account for 
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transcript isoforms), we reached an average transcriptome 
coverage of 23x for each replicate. An RPKM (Reads per 
Kilobase per Million reads) cutoff value of 0.1 was set to 
declare a locus expressed, resulting in 21,530 loci above the 
expression threshold. We also performed GO functional 
analysis to identify the biological processes that were most 
represented by the expressed genes. The GOSLIM terms 
for biological process assignment to the total expressed 
genes showed that cellular protein modification, carbohy- 
drate metabolic, transmembrane transport, translation and 
lipid metabolic processes were the most-represented bio- 
logical categories (Additional file 2: Figure S2). 

Differentially expressed genes and overrepresented 
categories 

To identify differentially expressed genes (DEGs), we 
used the R package DESeq to evaluate the significance of 
differences in expression and to control for false discov- 
ery rate (FDR). Pearson correlation coefficients among 
biological replicates were subjected to the same treat- 
ment and were all above 0.95. The FDR threshold was 
set to 0.1 and gene dispersion values were calculated by 
fitting a curve according to the DESeq "fit-only" mode. 
Using these criteria, we identified 712 DEGs (Additional 
file 3: Table S3), including 458 up-regulated genes and 
254 down-regulated genes. We also used MA plot ana- 
lysis to examine the magnitude distribution of the sig- 
nificantly regulated genes comparing the expression 
level to the log-transformed fold-change between MYC 
and CONT (Additional file 2: Figure S3). Among the 
DEGs, the 5 most highly-expressed genes (regardless of 
the growth condition) were a fruit-specific protein, Lino- 
leate 9S-lipoxygenase B, the acid beta-fructofuranosidase 



and two heat shock proteins 90 and 70 (Table 2). These 
genes affect important fruit processes like ripening, stress 
responses, and lipid and carbohydrate biosynthesis. 

GO annotation of differentially expressed genes 

We used GO enrichment analysis to identify the major 
gene groups affected by mycorrhization. As in all RNA- 
Seq experiments, the variable length of transcripts may 
produce biases in the data [23]; therefore, we analyzed 
GO enrichment via the Goseq R package, which was 
specifically designed to limit bias in the data [24]. We 
found 81 GO terms that were over-represented (FDR < 
0.1) in response to mycorrhization (Figure 2). The bio- 
logical processes corresponding to the lowest p-values, 
in decreasing order of importance, were "photosynthesis" 
(p = 2.66E-43), "photosynthesis light harvesting" (p = 
4.48E-24), "protein-chromophore linkage" (p = 2.45E-17), 
"oxidation-reduction process" (p = 1.74E-13). Other re- 
levant over-represented categories included: "response 
to biotic stimulus", "transmembrane transport" and 
"carbohydrate metabolic process" (Figure 2). 

Gene regulation In tomato fruit of mycorrhizal versus 
fertilized plants 

Examination of GO terms suggested that a large part of 
the modulated tomato transcriptome is devoted to con- 
trol of the main metabolic processes (i.e. photosynthesis 
and photorespiration) and to ripening (from cell-wall to 
secondary metabolites). For example, we identified 78 
genes involved in oxidation-reduction processes, 54 in 
metabolic processes, 49 in photosynthesis, 43 with a 
catalytic activity, 30 involved in transport process, 24 
transcription factors, 20 with a role in carbohydrate 



Table 1 Read number and alignment summary 


Sample 


Total number of filtered 
reads 


Total 
alignments 


Total unique aligned 
reads 


Total aligned 
bases 


Uniquely aligned mapping 
reads 


C0NT_1 


14,660,414 


13,136,222 


12,790,191 


1,038,315,933 


12,552493 


C0NT_2 


16,303,483 


14,602,025 


14,229,313 


1,141,477,858 


13,973,208 


MYC_1 


12,053452 


1 0,748,794 


10,510,262 


860,932,723 


10,336,271 


MYC_2 


11,715434 


10,222,971 


9,939,712 


814,018,186 


9,740416 
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Table 2 List of the five most strongly expressed genes 



Gene ID 


ITAG2.3_descriptions 


Description SL2.40 


IVIean read number 


Logjratio MYC vs CONTR 


Solyc07g049140.2.1 


Metallocarboxypeptidase inhibitor 


Fruit-specific protein 


200,347 


-1,097 


Solyc01g0991 90.2.1 


Lipoxygenase 


Linoieate 9S-lipoxygenase B 


149,420 


-1,178 


Solyc03g08391 0.2.1 


Acid beta-fructofuranosidase 


Acid beta-fructofuranosidase 


59,651 


-1,635 


Solyc07g0541 70.2.1 


Pectinesterase 


Pectinesterase 1 


55,636 


-1,154 


Solyc06g036290.2.1 


Heat shock protein 90 




29,385 


1,066 


Solyc04g011 440.2.1 


Heat shock protein 70 isoform 3 


heat shock protein 


25,164 


0,826 



The listed genes represent the five most expressed genes among the DEGs. The selection was based on the mean number of reads in MYC and CONT fruits. 



metabolism, 19 stress-responsive genes, 11 hormone- 
related genes and 51 unknown proteins corresponding 
to 10.95%, 7.58%, 6.88%, 6%, 4.17%, 3.37%, 2.8%, 2.66%, 
1.52%, and 7.12% of the regulated gene core set, respect- 
ively. In the following sections, we focus on a subset of 
these genes. 

Photosynthesis- and photorespiration-reiated processes 

We identified a substantial number of genes involved in 
different steps of photosynthesis as up-regulated in 
MYC fruit (Table 3). Some of these up-regulated genes 
encode proteins of photosystem II and photosystem I, 
and others are involved in light-dependent reactions 
(e.g. cytochrome b6f plastocyanins, ATP synthases and 
NADP reductases). Widespread up-regulation was also 
observed for genes related to the light-independent reac- 
tions or the Calvin cycle: two Ribulose bisphosphate car- 
boxylases, one phosphofructokinase family protein, four 
fructose-l-6-bisphosphatases class 1 and three fructose- 
bisphosphate aldolases. 

We also observed up-regulation of genes involved in 
chlorophyll metabolism. Six genes involved in chloro- 
phyll biosynthesis and one gene involved in chlorophyll 
degradation (encoding Senescence-inducible chloroplast 
stay-green protein 2) were also found to be up- 
regulated. Genes involved in chlorophyll biosynthesis 
encoded proteins including: magnesium chelatase H 



subunit, magnesium-protoporphyrin IX methyltransfer- 
ase, magnesium-protoporphyrin IX monomethyl ester, 
light dependent NADH:protochlorophyllide oxidoreduc- 
tase 1, chlorophyllide a oxygenase and geranylgeranyl re- 
ductase. In this context, it was of particular interest to 
observe the up-regulation of the Golden 2-like transcrip- 
tion factor (SIGLK2) (log2 ratio of 1.5) (Additional file 3: 
Table S3). This gene belongs to the GARP family of 
MYB transcription factors, and determines chlorophyll 
accumulation and distribution in developing fruit, thus 
influencing the entire photosynthetic process [25]. For 
chlorophyll degradation, we noted a slight up-regulation 
of the senescence-inducible chloroplast stay-green protein 
1 (SlSGRl) in MYC fruit (log2 ratio of 1). Stay-green genes 
encode members of a family of chloroplast proteins, in- 
volved in the brealcdown of chlorophyll-apoprotein com- 
plexes [26]. 

Genes involved in photorespiration (encoding serine 
hydroxymethyltransferase, hydroxypyruvate reductase, 
serine-glyoxylate aminotransferase, glycine decarboxyl- 
ase and glycolate oxidase were also up-regulated 
(Additional file 3: Table S3) in MYC fruit, suggesting that 
the photorespiratory activity in the fruit was probably 
higher under the mycorrhizal conditions. 

Taken as a whole, our transcript profiling data are 
consistent with numerous studies reporting the expres- 
sion of genes involved in light and dark photosynthetic 




letOa tet12 1et2t 1et30 1et39 le,03 1e,07 le+11 letlS 1et19 1et03 1et09 1et15 1e+21 1e+27 

1/(G0 term enrichment p value) 1/(GO term enrichment p value) 1/(GO term enrichment p value) 



Figure 2 GO enrichment in DEGs. The 81 enriched GO terms (y axis labels; FDR threshold < 0.1) associated to DEG, as assessed compared 
with all expressed genes (RPKM > 0.1), are plotted along the reciprocal of Enrichment p value (x axis) as calculated with the Goseq R package. 
Panel A: Biological Process; panel B: Molecular Function; panel C: Cell Compartment. 
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Table 3 Mycorrhiza-regulated genes involved in photosynthesis and related metabolisms in turning tomato fruit 



ITAG2.3_descriptions 



Description SL2.40 



Gene ID 



Log2 ratio IWYC vs 
CONTR 



Photosystem I 

Chlorophyll a/b binding protein 

Chlorophyll a-b binding protein 3C-like 

Chlorophyll a-b binding protein 37, 
chloroplastic 

Chlorophyll a/b binding protein 

Chlorophyll a/b binding protein 

Photosystem I reaction center subunit VI-1, 
chloroplastic 

Chlorophyll a-b binding protein 13, 
chloroplastic 

Chlorophyll a/b binding protein 

Photosystem I reaction center subunit 

Chlorophyll a/b binding protein 

Chlorophyll a-b binding protein 5A, 
chloroplastic 

Chlorophyll a-b binding protein 4, 
chloroplastic 

Photosystem I reaction center subunit V 

Chlorophyll a-b binding protein 8, 
chloroplastic 

Photosystem I reaction center subunit XI 

Photosystem I reaction center subunit III 

Photosystem I reaction center subunit VI-1, 
chloroplastic 

Photosystem I reaction center subunit IV A 
Photosystem I reaction center subunit XI 
Photosystem I reaction center subunit III 
Photosystem I reaction center subunit II 
Photosystem I reaction center subunit IV A 
Photosystem II 

Chloroplast photosystem ll-associated protein 

Photosystem II reaction center psb28 protein 

Photosystem II reaction center W protein 

Chlorophyll a-b binding protein, chloroplastic 

Chlorophyll a-b binding protein, chloroplastic 

Chloroplast photosystem II subunit X 
(Fragment) 

Oxygen-evolving enhancer protein 2, 
chloroplastic 

Oxygen-evolving enhancer protein 3 

Photosystem II reaction center W protein 

Oxygen-evolving enhancer protein 1 of 
photosystem II 

Oxygen-evolving enhancer protein 1 of 
photosystem II 



Chlorophyll a-b binding protein ID 
Chlorophyll a-b binding protein 3C 
Chlorophyll a-b binding protein 5 

Chlorophyll a-b binding protein IB 
Chlorophyll a-b binding protein IB 

Chlorophyll a-b binding protein 13 
Chlorophyll a-b binding protein IB 

Chlorophyll a-b binding protein IB 
Chlorophyll a-b binding protein 6A 

Chlorophyll a-b binding protein 4 
Chlorophyll a-b binding protein 8 



Photosystem I reaction center subunit II 



Photosystem II 22 kDa protein 



Chlorophyll a-b binding protein CP24 10A 
Chlorophyll a-b binding protein CP24 lOB 



Oxygen-evolving enhancer protein 2 



Solyc02g07 1000. 1.1 
Solyc03g005780.1.1 
Solycl2g0061 40.1.1 

Solyc02g070950.1.1 
Solyc02g071 01 0.1.1 
Solyc06g066640.2.1 

Solyc12g01 1450.1.1 

Solyc02g071 030.1.1 
Solyc08g01 3670.2.1 
Solyc02g070940.1.1 
Solyc05g056070.2.1 

Solyc07g047850.2.1 

Solyc07g0661 50.1.1 
Solyc10g007690.2.1 

Solyc06g082950.2.1 
Solyc02g069460.2.1 
Solyc03g1 20640.2.1 

Solyc06g083680.2.1 
Solyc06g082940.2.1 
Solyc02g069450.2.1 
Solyc06g054260.1.1 
Solyc09g063 130.2.1 

Solyc06g060340.2.1 
Solyc09g064500.2.1 
Solyc06g084050.2.1 
Solyc01g105030.2.1 
SolycOlgl 05050.2.1 
Solyc05g025600.1.1 

Solyc07g044860.2.1 



Photosystem II oxygen-evolving complex protein 3 Solyc02g079950.2.1 

Solyc09g06591 0.1.1 
Solyc02g090030.2.1 



Oxygen-evolving enhancer protein 1 



Solyc02g065400.2.1 



Inf 
5,75 

4,11 

3,35 
3,22 
2,92 

2,76 

2,68 
2,63 
2,61 
2,61 

2,54 

2,35 
2,23 

1,61 
1,31 
1,23 

1,10 
2,48 
2,41 
2,02 
1,73 

5,06 
2,53 
2,37 
2,19 
2,07 
1,97 

1,58 

1,16 
0,84 
1,51 

1,24 
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Table 3 Mycorrhiza-regulated genes involved In photosynthesis and related metabolisms in turning tomato fruit 

(Continued) 



Chlorophyll biosynthetic process 

Geranylgeranyl reductase 

Magnesium-protoporphyrin IX monomethyl 
ester 

Magnesium chelatase H subunit 

Magnesium-protoporpliyrin ix 
metliyltransferase 

Light dependent NADH:protocliloropiiyllide 
oxidoreductase 1 

CInloropliyllide a oxygenase 

Calvin cycle 

Ribulose bispliospiiate carboxylase small chain 
Ribulose bisphosphate carboxylase small chain 
Calvin, others 

Phosphofructokinase family protein 

Fructose-1 5-bisphosphatase class 1 

Fructose-1 5-bisphosphatase class 1 

Fructose-1 5-bisphosphatase class 1 

Fructose-1 5-bisphosphatase class 1 

Fructose-bisphosphate aldolase 

Fructose-bisphosphate aldolase 

Fructose-bisphosphate aldolase 

Chlorophyll breakdown 

Senescence-inducible chloroplast stay-green 
protein 2 

Cytochrome b6f complex and plastocyanin 

Plastoquinol-plastocyanin reductase 
Plastocyanin 

Cytochrome b5-f complex iron-sulfur subunit 
ATP synthases 

ATP synthase gamma chain 
ATP synthase delta subunit 
ATP synthase F1 delta subunit 
ATP synthase subunit b&apos 
NADP reductase 
Ferredoxin-NADP reductase 
Ferredoxin-NADP reductase 



Geranylgeranyl reductase 



Light dependent NADH:protochlorophyllide 
oxidoreductase 1 



Solyc03gn5980.1.1 
Solyc10g077040.1.1 

Solyc04g01 5750.2.1 
Solyc03g1 18240.2.1 

Solyc12g013710.1.1 

Solyc05g0503 10.2.1 



Ribulose bisphosphate carboxylase small chain 1 Solyc02g053 150.2.1 
Ribulose bisphosphate carboxylase small chain 2A Solyc03g034220.2.1 



Solyc07g045 150.2.1 
Solyc09g01 1810.2.1 
Solyc04g071 340.2.1 
Solyc10g085730.1.1 
Solyc05g052500.2.1 
Solyc02g052340.2.1 
SolycOlgl 10350.2.1 
Solyc02g084440.2.1 



Chloroplast sedoheptulose-1 ,7-bisphosphatase 



Green flesh protein; Senescence-inducible chloroplast Solyc08g080090.2.1 
stay-green protein 1 



Plastocyanin 



Solyc01g109040.2.1 
Solyc04g08201 0.1.1 
Solyc12g005530.1.1 

Solyc02g080540.1.1 
Solyc12g055830.1.1 
Solyc05g050500.1.1 
Solyc05g055000.1.1 

Solyc02g0521 30.2.1 
Solyc02g08381 0.2.1 



2,20 
2,10 

1,11 
1,09 

2,43 

1,92 

2,52 
2,25 

1,53 
3,15 
1,45 
1,35 
1,08 
3,23 
1,83 
1,17 

1,02 



2,59 
1,55 
1,39 

2,54 
1,43 
1,27 
0,85 

1,20 
0,98 



reactions, and photorespiration at turning and later rip- 
ening stages of the tomato fruit [13,27-29]. Biochemical 
data showed a decrease in enzymatic activity of glycolate 
oxidase and hydroxypyruvate reductase from breaker to 
red ripe stages [27] while a time-course analysis of 
genome-wide gene expression [13] showed that the large 
majority of photosynthetic and respiratory genes identi- 
fied in our experiment are down-regulated during the 
transition from green to red ripe fruit. Given that the 



nutritional status of the fruits from MYC and CONT 
plants was comparable, and all these genes were up- 
regulated in the MYC fruit, we hypothesize that mycorrhi- 
zation may help extend the duration of the photosynthetic 
and photorespiratory activities in tomato fruit. 

Carbohydrate and cell wall metabolism 

Tomato fruit ripening involves accumulation of sugars and 
breakdown of cell wall carbohydrates, resulting in fruit 
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softening. We identified twenty differentially regulated 
genes belonging to the carbohydrate metabolic process 
group, eight that are up-regulated in MYC fruit and twelve 
that are down-regulated (Table 4). Almost all genes in- 
volved in cell-wall carbohydrate modification were down- 
regulated (Figure 3) and many of them (glycoside hydrolase 
family 28 protein, beta-l-3-glucanase, beta-galactosidase, 
endo-l-4-beta-xylanase and polygalacturonase 2) are known 
to be involved in the degradation of cell-wall integrity, indu- 
cing fruit softening during ripening [30-33]. In the same cat- 
egory, the 5-xyloglucan endotransglucosylase/hydrolases 
iXTH) which encode enzymes with xyloglucan endotrans- 
glucosylase (XET) or endohydrolase (XEH) activities were 
also down-regulated. The role of this gene family during 
tomato ripening has been extensively investigated [34-36] 
and XTH genes show different behaviors during ripening. In 
particular, XTHl, XTH7, XTH12 and XTH6 had their high- 
est expression in the mature green stage and then decreased 
in the later stages. 

Our data also show the down-regulation of other 
genes encoding cell wall associated proteins such as 



pectinesterases and expansins (Table 4). Expansins are 
involved in expansion of plant cells in growing tissues 
but were also identified in non-growing tissues [37]. One 
of the expansins identified in our experiment is LeExpl 
(Solyc06g051800.2.1) which was found to be specifically 
expressed in ripening fruit, with expression that was very 
weak at the green stage and increased with ripening 
[13]. In addition, LeExpl seemed to be directly regulated 
by ethylene at the transcriptional level [37]. The second 
class of cell-wall modifying proteins identified was the 
pectinesterases. Pectinesterases demethylate pectins and 
make them more susceptible to polygalacturonase action 
[38], but they may also reinforce cell-to-cell adhesion 
[39]. The down-regulated pectinesterase gene identified 
in our experiment {PME 1.9) (Solyc07g064170.2.1) was 
one of the five most highly-expressed genes among the 
DEGs (Table 2). 

All together, the results related to the regulation of 
genes encoding cell wall modifying enzymes show that, 
with few exceptions, the genes involved in cell wall hy- 
drolysis are repressed in the MYC fruit, suggesting that 



Table 4 Mycorrhiza-regulated genes involved in carbohydrate and cell-wall metabolism In turning tomato fruit 



Gene ID 



ITAG2.3_descriptions 



Description SL2.40 



Log2ratio IVIYC vs 
CONTR 



Carbohydrate metabolism 

Solyc10g083290.1 .1 Beta-fructofuranosidase insoluble isoenzyme 2 



Solyc09g011 81 0.2.1 
Solyc02g082920.2.1 
Solyc11g007990.1.1 
Solyc04g071 340.2.1 
Solycl0g086730.1.1 
Solyc05g052600.2.1 
SolycOSgOl 381 0.2.1 
Solyc03g1 13030.2.1 
Solyc08g082 170.2.1 
Solyc04g016470.2.1 
Solyc02g084720.2.1 
Solyc03g08391 0.2.1 
Solyc07g056000.2.1 
Solyc03g093 130.2.1 
SolycOl gill 340.2.1 
Solyc01g099630.2.1 
Solyc03g093 120.2.1 
Solyc05g005 170.21 
Other cell wall genes 
Solyc07g0641 70.2.1 
Solyc03gl 23630.21 
SolycOl g1 12000.21 
Solyc06g05 1800.2.1 



Fructose-1 6-bisphosphatase class 1 
Endochitinase (Chitinase) 
Malate dehydrogenase 
Fructose-1 6-bisphosphatase class 1 
Fructose-1 6-bisphosphatase class 1 
Fructose-1 6-bisphosphatase class 1 
Glycosyl hydrolase family 5 protein/cellulase 
Aldse 1-epimerase-like protein 
Glycoside hydrolase family 28 protein 
Beta-1 3-glucanase 
Beta-galactosidase 
Acid beta-fructofuranosidase 
Xyloglucan endotransglucosylase/hydrolase 7 
Xyloglucan endotransglucosylase/hydrolase 9 

Endo-1 4-beta-xylanase 
Xyloglucan endotransglucosylase/hydrolase 5 
Xyloglucan endotransglucosylase/hydrolase 9 
Polygalacturonase 2 

Pectinesterase 
Pectinesterase 
Expansin-like protein 
Expansin 



Cell-wall invertase 

Acidic 26 kDa endochitinase 
Malate dehydrogenase 

Chloroplast sedoheptulose-1 ,7-bisphosphatase 



Glucan endo-1, 3-beta-D-glucosidase 
Beta-galactosidase 
Acid beta-fructofuranosidase 
Xyloglycan endo-transglycosylase 
Xyloglucan endotransglucosylase-hydrolase XTH3 

Probable xyloglucan endotransglucosylase/hydrolase 1 



Pectinesterase 1 
Pectinesterase/pectinesterase inhibitor U1 
Expansin-like protein 
Expansin 1 protein 



4,35 
3,16 
1,76 
1,66 
1,46 
1,35 
1,08 
1,06 
-0,82 
-0,83 
-1,31 
-1,40 
-1,63 
-1,70 
-2,28 
-2,33 
-2,37 
-2,49 
-5,01 

-1,15 
-1,34 
-1,53 

-1,62 
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Figure 3 Clustering and heatmap of expression values for DEGs in the cell wall structure-related gene class. Light and dark blue indicate 
Iniglier and lower expression values, respectively. Gene descriptions are as reported in Tomato Genome Consortium (2012). The heatmap was 
generated starting from DESeq normalized counts and further processed with custom scripts based on heatmap.2 function as available in the 
'gplots' Bioconductor package and transformation function VST (DESEq package). 



the molecular mechanisms that lead to softening are de- 
layed in the MYC fruit. 

The MYC fruit also showed down-regulation of an 
acid beta-fructofuranosidase (Solyc03g083910.2.1; log2 
ratio: -1.635) coupled with a higher up-regulation of a 
cell-wall invertase (SolyclOg083290.1.1; logs ratio: 4.35) 
(Table 4). BLAST searches of these two sequences showed 
that they are, respectively, the tomato vacuolar invertase 
TIVl and the cell -wall invertase {LIN6). These vacuolar 
and cell wall acid invertases, cleave the sucrose transported 
from source organs (leaves) into hexoses (glucose and fruc- 
tose). Heterotrophic organs like fruit then use the hexoses 
as a direct source of carbon and energy [40]. The supply of 
sucrose to fruit and other sink tissues can limit their 
growth [41]. In carbohydrate supply to sink organs, sieve 
elements release sucrose into the apoplast where it under- 
goes an irreversible hydrolysis catalyzed by extracellular in- 
vertase and sucrose synthase [42]. The resulting hexoses 
are then transferred by hexose transporter(s) into sink cells. 
In tomato, extracellular invertases are encoded by LIN5, 
LIN6, LIN7 and LIN8, which have highly tissue-specific 
expression, regulated by internal and external signals [42]. 
TIVl is expressed in red tomato fruit and LIN6 is 
highly expressed in seedling roots, flowers and 



tumors, but not in fruit [42]. Our RNA-Seq data show 
very low expression of LIN6 in CONT fruits (1.35 nor- 
malized read counts) but a much higher expression in 
MYC fruits (27.75 normalized read counts). Together, 
these data suggest that sucrose is actively mobilized in 
the MYC fruit, providing ready-to-use sugars for fruit 
development. 

Transport processes 

We identified 30 DEGs as involved in transport pro- 
cesses (Figure 4), and some of them are related to the 
DEGs involved in primary and carbohydrate metabol- 
ism. Among them, one of the most up-regulated genes 
in MYC fruit was a hexose transporter 2 {LeHT2) 
(Solyc09g075820.2.1; logj ratio: 3.24). This gene shows 
the highest sequence similarity to a glucose/H + and 
low-affinity fructose/H + symporter that is expressed in 
tomato fruit storage parenchyma cells and strongly in- 
duced at fruit maturity [43,44]. Hexose transporters 
transfer the hexoses hydrolyzed by cell wall invertases 
into sink cells. Therefore, they play an important role in 
sugar accumulation by fruit storage cells [45]. Indeed, a 
significant reduction (80-90%) in expression levels of 
three hexose transporter genes, LeHTl, LeHT2, LeHT3, 
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Figure 4 Clustering and heatmap of expression values for DEGs in the transmembrane transport-related gene class. Light and dark blue 
indicate higher and lower expression values, respectively. Descriptions as obtained from BLAST2G0 annotation are shown. The heatmap was 
generated starting from DESeq normalized counts and further processed with custom scripts based on heatmap.2 function as available in the 
'gplots' Bioconductor package and transformation function VST (DESEq package). 



by RNAi-mediated knockdown caused a 55% decrease 
in fruit hexose accumulation [44] . 

A tomato hexose transporter has been described as 
up-regulated after infection with Tomato yellow leaf curl 
virus and as co-regulated with a tomato lipocalin [46]. 
Lipocalins transport small hydrophobic molecules, and 
play important roles in various biological processes such 
as development and stress responses. Interestingly, a to- 
mato lipocalin (Solyc07g005210.2.1) was activated in the 
MYC fruits, suggesting a network of interactions be- 
tween transporters and responses to biotic stresses (see 
below). 

We observed the up-regulation of one plasma mem- 
brane intrinsic protein and three tonoplast intrinsic pro- 
tein aquaporins. Aquaporins function in water transport, 
and play pivotal roles in regulating water movements in 
fruit, which are crucial for fruit expansion and ripening 
[47] . Aquaporins are highly responsive to mycorrhization 
in legumes [48-52], although their exact role in the sym- 
biosis is not clear [53] and tomato roots may also show 
a downregulation of aquaporins [6] . 

A urea transporter, characterized by a high up-regulation 
in the MYC fruit, showed 71% amino acid sequence simi- 
larity with an Arabldopsis thaliana sodium/urea cotran- 
sporter (DUR3). Interestingly, an aUantoinase gene was 
found to be up-regulated in MYC fruit of the Micro-Tom 



variety, suggesting that the metabolism of urea precursors 
is also activated [8]. 

In accordance with the prominent function of trans- 
port in mycorrhizal roots, the membrane transporters 
functional category was found to be one of the most- 
regulated categories in roots colonized by AM fungi 
[50,51,54]. 

Our experiments identified two up-regulated genes in- 
volved in phosphate transport. The inorganic phosphate 
transporter (Solyc09g066410.1) has homology to both 
LePTl and LePT3 phosphate transporters. Unlike LePTl, 
which is constitutively expressed [4], Solyc09g066410.1 
seems to be induced by mycorrhization as it showed no 
basal expression in the iruit of non-inoculated plants. The 
second identified phosphate transporter (Solyc05g013510.2) 
shows the highest similarity with a chloroplastic inorganic 
phosphate transporter from S. tuberosum; this gene was 
highly expressed in young, rapidly growing tissues and 
may have a role in cell wall metabolism [55]. These data 
provide a good mechanistic basis for the observed fruit P 
content, which reaches comparable values in the two 
MYC and CONT conditions. This suggests that the 
mycorrhizal- dependent P uptake is as efficient as the P 
uptake from fertilization. Similarly, two sulphate-related 
transporters (High affinity sulphate transporter, see 
Additional file 3: Table S3) were more strongly 
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expressed in the fruits from mycorrhizal plants, consist- 
ent with the comparable S content. 

A solute carrier protein, probably acting as a nitrate 
transporter, was found to be down-regulated in our 
dataset. Nitrate transporters were found to be either 
up-or down-regulated in Medicago mycorrhizal roots, 
depending on the specific gene considered [48,51,56]. 
In tomato, colonization by F. mosseae specifically 
enhanced the expression of the nitrate transporter 
LeNRT2;3 without affecting the expression of other 
isoforms [57]. 

In conclusion, in addition to the expression of genes 
that function in fruit development, the presence of 
the AM fungus activates genes involved in nutrient 
and water transport activities, which are strictly re- 
lated to the known functional markers of mycorrhizal 
roots. 

Regulation of hormonal pathway genes 

Tomato fruit development and ripening is stringently 
regulated by hormone dynamics [22]. In roots, the hor- 
monal balance modulates mycorrhization and vice- 
versa [7,58-60]. In addition, tomato genes involved in 
hormonal metabolism show a systemic response to the 
AM fungus [6]. We therefore expected hormone path- 
way genes to be differentially regulated in fruit; indeed, 
11 genes were found to be involved in gibberellin, absci- 
sic acid (ABA) and ethylene metabolism (Additional file 
3: Table S3). A Gibberellin 2-oxidase 1 involved in the 
degradation of the active gibberellic acid (GA), and two 
Gibberellin regulated proteins, GASTl (gibberellic acid 
stimulated transcript) and RSI-1 (Root system inducible- 
1), were found to be down-regulated in the fruit of MYC 
plants. Both GASTl and RSI-1 belong to the GASA (GA- 
stimulated Arabidopsis) gene family, which includes many 
hormone-responsive genes involved in pathogen responses 
as well as in different aspects of plant development [61]. 
Two genes involved in ABA metabolism were found to be 
up-regulated in MYC fruit: the first one is zeaxanthin 
epoxidase, which catalyzes the conversion of zeaxanthin 
to violaxanthin [62], an ABA precursor [63]. The second 
gene is the ABA 8 -hydroxylase CYP707A2, which plays a 
role in ABA catabolism and is activated by the presence of 
ABA [64]. 

Six ethylene responsive factors were down-regulated in 
the MYC fruit (Additional file 3: Table S3): Ethylene re- 
sponse factor 1 {ERFl), a probable ethylene responsive 
transcription factor 2b, the C-repeat-binding factor 3 
(CBFS), a probable ethylene-responsive transcription fac- 
tor 1, the ethylene-responsive transcription factor 10 
and the ethylene-responsive late embryogenesis-like pro- 
tein. The ERFs are secondary ethylene activated factors, 
which increase plant tolerance to stress by modulating 
expression of downstream stress-responsive genes [65]. 



ERFl is a key element in the response to different necro- 
trophic pathogens [66]. Interestingly, in Medicago the 
transcription factor ERF 19 is down-regulated by the per- 
ception of an effector molecule secreted by an AM fun- 
gus [67]. However its most closely related homologue in 
tomato does not show differential expression during 
fruit ripening. 

Overall, these data demonstrate the dynamic expres- 
sion of genes involved in hormone metabolism: on the 
one hand, they mirror the transcriptional profile of 
fruit ripening markers (ethylene-responsive genes) and 
related pathways; on the other, some of them mirror 
the known role of hormone networks in AM estab- 
lishment. Indeed, a positive correlation between ABA 
and mycorrhization in roots was reported [58] and 
ABA was found to be necessary for complete and fully 
functional arbuscule development and for successful 
colonization of the plant root [58]. By contrast, nu- 
merous studies have shown a reciprocal negative rela- 
tionship between ethylene and mycorrhization [7,59]. 
However, the role of hormone networks in AM estab- 
lishment remains far from fully understood [7,60]. 

Terpenoid metabolism 

Terpenoids include phytohormones (abscisic acid, gib- 
berellins, cytokinins, and brassinosteroids), photosyn- 
thetic pigments (phytol and carotenoids), and electron 
carriers (ubiquinone) [68]. However, most terpenoids are 
specialized secondary metabolites that function as defen- 
sive compounds accumulating upon pathogen infection 
(i.e. phytoalexins) [69] or attracting predators/parasitoids 
of the attacking insect [70]. In tomato, terpenoids are 
present in large quantities in the glandular trichomes on 
leaves, stems, young fruits, and flowers [71-74]. 

Three genes belonging to the terpene synthase [TPS) 
family, TPS31, TPS32 and TPS33, were induced in the to- 
mato fruit by mycorrhization. Although only a few genes 
were affected, these three genes were among the ten most 
up-regulated by mycorrhization (Additional file 3: Table S3) 
and TPS31 showed no basal expression in CONT tomato 
fruits. TPS31, TPS32 and TPS33 belong to the sesquiterpe- 
noid phytoalexin biosynthesis pathway (Additional file 2: 
Figure S4), leading to the production of the phytoalexin 
rishitin, well known for its antimicrobial properties [69]. 
Accumulation of phytoalexins has been suggested as a 
defense mechanism in mycorrhizal plant roots [75-77]. 
By contrast, genes involved in carotenoid synthesis did 
not show any change in their transcriptomic profiles 
(Additional file 2: Figure S5). These data were con- 
firmed by a biochemical analysis of these metabolites, 
which showed no statistically significant differences be- 
tween the two sets of fruits (Additional file 1: Table S4), 
with the exception of general reducing activity, which 
was higher in the MYC fruits. 
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Defense genes 

Plants activate a wide range of mechanisms to defend 
themselves against external attacks [78]. Such responses 
are under the control of a finely tuned hormonal network, 
and can spread from the infection site to a systemic level. 
In our study, four out of the five up-regulated genes were 
pathogenesis related proteins (PRs) three PRlOs, and one 
PRl. PRs specifically accumulate in plants in response to 
pathogen infection and act in development of systemic ac- 
quired resistance against further infections. PR regulation 
in response to mycorrhization appears to be complex; 
Lopez-Raez and colleagues [7] reported the induction of 
PRIO in tomato roots inoculated with Funneliformis mos- 
seae and Rhizophagus intraradices. By contrast, some PRs 
were down-regulated in plant shoots following mycorrhi- 
zation [6,49] and induced only after a pathogen attack, as 
a result of the mycorrhizal priming effect [79]. In our ex- 
periment, we observed that mycorrhization per se induced 
the up-regulation of PRs in tomato fruit (Additional file 3: 
Table S3), as already observed for PRl in rice leaves [80]. 
Another gene (an acidic endochitinase) up-regulated in 
the fruit from MYC plants can be grouped in this cat- 
egory, as chitinase induction results in the activation of a 
mycorrhiza-mediated defense response to worms in grape- 
vine roots [81] as well as in mycorrhizal rice leaves [80]. 

Finally, a gene coding for a plant defensin containing 
a gamma-thionin domain was found to be down- 
regulated. Defensins are part of the front line of the 
plant immune system, inhibiting the growth of bacteria 
and fungi. Strikingly, three defensins were highly in- 
duced in AM-colonized Medicago roots, and their 
expression was strongly up-regulated in arbuscule- 
containing cells [51]. 

A number of papers have reported expression of ho- 
mologues of these tomato genes in roots, shoot and 
leaves of mycorrhizal plants, opening the question 
whether AM fungi influence host immunity by modulat- 
ing defense responses [82]. Our data suggest that some 
of these responses also remain active at the systemic 
level. In conclusion, it seems that a set of defense- 
related genes -already identified in mycorrhizal plants - 
are differentially expressed in tomato fruit at the turning 
stage, suggesting that they could take part in the so- 
called priming reaction elicited by mycorrhizal fungi and 
are also maintained at systemic level in the fruit. 

Time course analysis and validation of selected DEGs 

To validate the RNA-Seq data and to check whether the 
effect of mycorrhization on gene expression maintained 
the same trend (up-or down-regulation) at different stages 
of ripening, we used quantitative real-time PGR (q-PCR) 
to examine the pattern of gene expression of nine selected 
DEGs at three fruit ripening stages: mature green, breaker 
and turning. We also tested a group of non-regulated 



genes involved in carotenoid metabolism: 1-deoxy-D-xylu- 
lose-5-phosphate synthase {DXS), 4-hydroxy- 3-methylbut- 
2-en-l-yl diphosphate reductase {HDR), phytoene synthase 
(PSY), carotenoid cleavage dioxygenase (CCD), and lipoxy- 
genase (LOX). The trend of up- or down-regulation ob- 
served in the RNA-Seq experiment was validated for aU 
the genes at the turning stage, and statistically significant 
changes were observed for 7 out of 9 genes in at least one 
time point (Figure 5; Additional file 2: Figure S6). None of 
the carotenoid metabolism genes turned out to be regu- 
lated by mycorrhization in q-PCR, according to the RNA- 
seq results (Additional file 2: Figure S5). These results 
verify the reproducibility and reliability of the RNA-Seq 
technique and suggest that the regulation of the selected 
genes by mycorrhization is not restricted to the turning 
stage at which the RNA-Seq was performed. Mycorrhiza- 
tion often exerted a more intense effect at other ripening 
stages, when a peak in gene expression was observed 
(Figure 5). The values were, in most cases, higher in RNA- 
Seq than in q-PCR experiments; this was mainly observed 
for genes with low fold-change levels. The time course q- 
PCR analysis allowed us to validate the changes in expres- 
sion for genes with particular importance. For instance, 
the statistically significant up-regulation of IPT in all the 
three ripening stages allowed us to hypothesize that this 
gene may be a mycorrhiza-induced transporter, like other 
tomato phosphate transporters already found to be acti- 
vated by mycorrhization in the roots [5,83]. Our analysis 
indicates that ZPr could be a good candidate to be used as 
a fruit marker for mycorrhization. 

Discussion 

Our genome-wide RNA-Seq analysis of tomato fruits 
from mycorrhizal plants demonstrated that the presence 
of the mycorrhizal fungus in the roots influences genes 
involved in some processes strictly related to the ripen- 
ing of the fleshy tomato fruit. Our strategy to compare 
mycorrhizal plants with fertilized plants allowed us to 
uncouple the nutritional effects from the presence of 
AM fungi, which improve plant nutrition [20]. In this 
way we avoided the stress conditions that are character- 
istic of control non-fertilized plants, as plants respond to 
nutrient deficiency with important changes in root and 
shoot architecture [84]. We cannot exclude that the dif- 
ferent fertilization treatment may have affected the tran- 
scriptomic profile; however, our data strongly suggest 
that, in fruits, AM symbiosis affects a core set of genes 
that have already been identified as sensitive to the pres- 
ence of AM fungi (i.e. nutrient transporters). Along the 
same line, AM symbiosis also affects transcripts involved 
in major fruit pathways including photosynthesis and 
photorespiration as well as amino acid and cell wall 
metabolism. Figure 6 shows an overview of the major 
expression changes that we interpreted as caused by 
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Figure 5 q-PCR time course analysis of gene expression in fruits of different stages. Relative gene expression in mature green (G), breaker 
(B) and turning (T) stage fruits, of nine selected differentially expressed genes (red line). The blue bar indicates the log2 fold change obtained at 
the turning stage in the RNAseq experiment. Identities of genes are as follows: Solyc01g008540.2.1 (Cinnamoyl CoA reductase-like protein); 
Solyc02g078400.2 (Allantoinase); Solyc02g076920.2.1 (Transcription factor MYC2); Solyc10g005080.2.1 (Late elongated hypocotyl and circadian 
clock associated-1-like); Solyc09g066410.1 .1 (Inorganic phosphate transporter); Solyc06g062540.2.1 (Putative acid phosphatase); Solyc03g08391 0.2.1 
(Acid beta-fructofuranosidase); Solyc01g080280.2.1 (Glutamine synthetase); Solyc10g055390.1 (Nodulin family protein). 



mycorrhization, as annotated by MapMan software [85]. 
The ft'uits of mycorrhizal plants show transcriptomic 
signatures characteristic of a climacteric fleshy fruit, and 
of the mycorrhizal status of its roots. Thus, even if bio- 
logical conclusions based on a transcriptome analysis re- 
main speculative, the obtained results, coupled with 
information on fruit nutrient contents, allow us to dis- 
cuss some hypotheses on the mechanisms by which AM 
fungi might influence tomato fruit. 

From roots to fruits: sharing similar transcription profiles 

The bulk of the data available on gene expression in 
mycorrhizal conditions relates to tomato roots [54] as 
the physical site of mycorrhizal interaction and function- 
ing. The fruit possesses a unique metabolism, leading to 
the production of specific metabolites, including pig- 
ments, flavor, and aroma compounds. It is thus particu- 
larly interesting, but at the same time tricky, to draw a 
parallel between the functional meaning of gene regula- 
tion in roots and what we observed for the first time in 
the fruit transcriptome of mycorrhizal tomato plants. 

Out of the 712 regulated genes in the MYC fruit, 5 genes 
have homologues that have already been reported as acti- 
vated in mycorrhizal plants; these include a TIP-aquaporin 
expressed in tomato leaves [86] a pathogenesis-related 
protein 10 [7] expressed in roots, and a probable root 



mycorrhiza-induced inorganic phosphate transporter [5]. 
Our RNA-Seq analysis also allowed us to validate the 
changes in expression of some genes identified in our pre- 
vious microarray experiment on Micro-Tom fruit as regu- 
lated by mycorrhization [8]. For example, we confirmed 
the up-regulation of the allantoinase gene and the down- 
regulation of the glycolate oxidase 2 gene [8]. 

A huge number of transcriptomic investigations have 
shown the extensive effects of AM colonization on host 
metabolism [51,54] and identified a core set of genes 
that respond to mycorrhization for all tested plants, 
fungi, and organs. Our global transcriptome profiling, 
which analyzed the tomato fruit transcriptome via RNA- 
Seq, has identified two gene categories where the sys- 
temic effect of the mycorrhization is of particular rele- 
vance: the transporters and the defense-related 
responses. 

The up-regulation of transporter genes illustrates the 
specific lifestyle of mycorrhizae; AM fungi support the 
plant, improving its mineral nutrition and receiving fixed 
carbon. The regulation of transporter-coding genes in 
organs other than the roots strongly suggests that this 
function is also highly responsive to mycorrhization at a 
systemic level. This leads to the major functional result 
of our analysis: the transcriptomic signature and the 
mineral contents of the fruits demonstrate that the AM 
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Figure 6 Overview of tlie effect of mycorrhlzation-lnduced expression cfianges on metabolism of tomato fruits. MapMan software 
(Metabolism_overview panel) was used to provide a snapshot of modulated genes over the main metabolic pathways. DEGs were binned to 
MapMan functional categories and Log2 fold changes values are represented. Up-regulated and down-regulated transcripts are shown in red and 
blue, respectively. 



fungus improves plant and fruit nutrition as efficiently 
as heavy fertilization. It would be of interest to under- 
stand whether in field conditions, mycorrhization could 
confer the same advantages as chemical fertilization, in 
terms of fruit quality and/or energy consumption. 

Similarly, several studies have examined the effects of 
mycorrhization on plant defense mechanisms. Recent re- 
search suggests that mycorrhizal establishment stimu- 
lates priming of the plant's innate immune system rather 
than the induction of specific defense mechanisms 
[82,87]. Mycorrhizal plants seem to react more promptly 
to pathogen attack, eventually showing enhanced resist- 
ance [80,88]. Our results on the activation of PR pro- 
teins, homologs of others already reported as activated 
in AM roots [7], show that in the MYC fruit, more com- 
plex regulation takes place, as some of these genes were 
down-regulated in shoots [6,49] while others were up- 
regulated [80]. 



The AM fungus in the root affects the transcriptomic 
profile of genes involved in fruit photosynthesis, 
photorespiration and sink strength 

Developing tomato fruits conduct photosynthesis and 
contribute up to 20% of the fruit's photosynthate, as 
shown by their expression of proteins involved in light- 
harvesting, electron transfer and CO2 fixation [89]. Dur- 
ing ripening, tomato fruits shift from partially photo- 
synthetic to truly heterotrophic metabolism, as also 
mirrored by a decrease of many photosynthetic proteins 
together with those involved in chlorophyll biosynthesis 
[29]. However, notwithstanding the wealth of data, the 
role of fruit photosynthesis remains not fully understood 
[90]. Our transcriptomic analysis adds a novel piece of 
information to this emerging understanding, showing 
that mycorrhization might affect fruit photosynthe- 
sis and photorespiration. This transcriptional influence 
seems to be independent of the nutritional level, since 
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mycorrhizal and control fertilized plants had similar 
mineral and nutrient contents in their fruits (Additional 
file 1: Table SI). 

Here we found that tomato fruits of mycorrhizal 
plants sampled at the turning stage present an important 
up-regulation of photosynthetic genes, which are usually 
activated in immature green fruits [91]. Also, in MYC 
fruits photorespiration may be activated via the in- 
creased expression of a core set of genes (Additional 
file 2: Figure S2). AM fungi stimulate photosynthesis in 
leaves of different plants [92,93]. Interestingly, strigolac- 
tones, which act as positive regulators of AM symbiosis, 
also up-regulate light-harvesting genes in the same to- 
mato organs [94]. It has been proposed that mycorrhi- 
zation causes an improvement of plant water status, 
increasing hydraulic conductance and affecting stomata 
and leaf expansion [95,96]. Other investigations sug- 
gested that the AM fungi, which depend on their host 
for carbon, increase root sink strength and enhance leaf 
photosynthetic activity to compensate for the 'cost' im- 
posed upon the plant C economy by the fungus [97,98]. 

On the basis of our RNA-Seq profiling, we hypothesize 
that similar events also take place in the fruit, where the 
increased sink strength is confirmed by the activation of 
other metabolic pathways. The up-regulation of LIN6 in 
our MYC fruit is particularly intriguing as the expression 
of this gene has been reported to be high in roots, flowers 
and tumors and was almost absent in fruit [42]. Indeed, in 
our experiment, LIN6 was very weakly expressed in the 
CONT fruit. 

It seems therefore that mycorrhization positively influ- 
ences the expression of genes of the metabolic pathways 
related to carbon assimilation and respiration. In this 
context, the up-regulation of the S1GLK2 transcription 
factor [25] could provide novel and unexpected insights. 
S1GLK2 is encoded by the uniform ripening locus (m) 
and determines the intensity and pattern of chlorophyll 
distribution in the unripe fruit. The dominant U allele 
leads to fruit with dark green shoulders at the stem end 
adjacent to the pedicel, and u/u fruit are uniformly light 
green and were therefore favored by the breeders. How- 
ever, in these fruits selection has led to suboptimal 
chloroplast development, and consequently, to decreased 
ripe fruit sugar and lycopene levels. It seems that in 
some way mycorrhization rescues the levels of SIGLK2 
mRNA to higher levels, also probably controlling the 
downstream photosynthesis-related genes. 

Could mycorrhization extend tomato shelf life? 

The ripening of climacteric fruits such as tomato is char- 
acterized by a large burst of respiration known as the 
'climacteric rise' at the onset of ripening. The climacteric 
rise is coupled with increased ethylene biosynthesis [39], 
which is thought to coordinate numerous biochemical 



and physiological processes necessary for normal ripen- 
ing such as softening of the cell wall, accumulation of 
pigment and increases in flavor and aroma [14]. Cell wall 
softening, an important process influencing fruit quality 
and principally shelf life, involves a coordinated decline 
in the activity of cell-wall biosynthetic enzymes and an 
increase in the activity of cell wall hydrolytic enzymes, 
like polygalacturonase [99]. 

Availability of the high-quality, annotated tomato gen- 
ome sequence [13] has allowed an extension of our under- 
standing of ripening, identifying more than 40 genes 
involved in fruit cell wall architecture. The regulation of 
cell wall related genes strictly relates phenotypic features, 
including fruit firmness, and shelf life. Our transcriptomic 
data consistently demonstrate that mycorrhization may 
lead to a down-regulation of cell wall hydrolytic enzymes. 
Similar results were already reported in a previous, more 
limited study [8] of the Micro-Tom variety. 

All these data, together with the up-regulation of the 
SIGLK2 transcription factor, the up-regulation of photo- 
synthesis-related genes, the increased chlorophyll synthe- 
sis and the down-regulation of ethylene responsive genes, 
indicate that mycorrhization may lead to fruits with a po- 
tentially longer shelf-life through the activation of multiple 
pathways. 

Conclusions 

The extensive transcriptomic analysis we performed, com- 
paring tomato fruits from mycorrhizal plants with fruits 
from fertilized plants, demonstrates that mycorrhization 
may have important systemic effects. By improving the 
nutritional status of the whole plant, mycorrhization leads 
to a mineral content in the fruits that is comparable to the 
mineral content of fruits from highly fertilized plants. By 
affecting the source-sink relationships of the whole plant, 
mycorrhization also has an impact on some of the features 
of fleshy fruits. This novel transcriptome survey thus 
shows biological intersections between roots and fruits, 
and will provide a platform for future, in-depth studies. 
The data we present here will allow study of the complex 
process of tomato fruit development under the usual culti- 
vation conditions, like those occurring in greenhouses 
where the plants are fertilized, the soil is not sterile, and 
the complex tomato microbiome, also including AM 
fungi, affect plant performance. 

Methods 

Plant materials and growth conditions 

Solarium lycopersicum cv. Moneymaker tomato plants 
were grown in the greenhouse from March through De- 
cember under natural light, which provided a photoperiod 
ranging from 9 h to 15 h, temperature from 16°C to 30°C 
and relative humidity from 55% to 80%. The seedlings 
were grown in pots with sterile quartz sand (80%) and 
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Sterile pumice (20%) for non-inoculated plants or with a 
mixture of sand (55%), sterile pumice (20%) and Funne- 
liformis mosseae (formerly Glomus mosseae) inoculum 
(25%) Gerd. & Trappe (BEG 12) purchased from Agraux- 
ine (Bretagne, France, http://www.agrauxine.fr) for inocu- 
lated plants. Forty-five plants divided into three groups of 
fifteen plants according to their growth conditions were 
considered: plants grown with F. mosseae inoculum 
(MYC) or without inoculum were watered, once a week, 
with a modified Long-Ashton solution containing a low 
phosphorus and nitrogen concentration (3.2 i^M P, 
1 mM N, 1.75 mM S, 1 mM K) [100]. The fertilized con- 
trol plants (CONT) were grown without inoculum and 
watered, once a week, with a commercial fertilizer solution 
(3.8 mM P, 12.9 mM N, 6.6 mM S, 3.8 mM K) (Asso di 
fiori, CIFO S.p.A). Phenological features were monitored 
during the experiments, as summarized in Additional file 1: 
Table S2. At the end of the experiment, roots from MYC 
and CONT plants were cut and fungal colonization was 
assessed according to the Trouvelot system [101] using 
MYCOCALC software. 

Illumina sequencing 

For the RNA-Seq experiment, fruits were harvested from 
the CONT and MYC plants when the fruits reached the 
turning stage. RNA was extracted using the 'pine tree- 
method' [102] with the addition of 1% PVPP to the ex- 
traction buffer. For each growth condition, we used two 
biological replicates, each containing the pooled RNA 
from three fruits. Ten micrograms of each RNA sample 
was sent to Fastens Life Science Co. (Geneva, Switzerland) 
where the libraries were produced and sequenced using 
the Illumina Genome Analyzer (Solexa). The four libraries 
were indexed and single-end multiplexed sequencing was 
performed using 100 bp length reads. 

Time course quantitative Real-Time PCR validation 

Quantitative Real-Time PCR, was used to measure the 
expression of nine genes shown to be differentially regu- 
lated by RNA-Seq, as well as of five carotenoid related 
genes. Three fruit ripening stages were considered: ma- 
ture green, breaker and turning. For each stage, three 
biological replicates were used. RNA extraction, primer 
design, quantitative real time PCR experiment and data 
analysis were carried out as described in Salvioli et al. 
[8]. The primer names and corresponding sequences are 
listed in Additional file 1: Table S5. 

Blolnformatic methods 
Mapping of Illumina reads 

Raw fastQ files were checked for contaminants and low- 
quality bases and contaminants were removed with the 
cutadapt software [103]. The spliced read mapper TopHat 
version 1.4.1 [104] was used to map reads to Tomato 



SL2.40 (ITAG2.3) obtained firom the plant repository. A 
minimum and maximum intron length of 40 and 50,000 
bp were used, respectively. Read counts were collected 
with HTSeq version 0.5.3 (http://www-huber.embl.de/users/ 
anders/HTSeq) in the single end and 'union' mode using So- 
larium lycopersicum MSU SL2.40-16 GTF file as obtained 
from the ensembl plant repository. 

DEG calling 

The DESeq Bioconductor package version 1.10.1 [105] 
was used to call differentially expressed genes. DESeq im- 
plements a model based on negative binomial distribution; 
this model was developed with special attention to coping 
with biological variance and was run under R release 
2.15.2. The cutoff for considering a gene expressed was set 
to 0.1 RPKM (Reads per Kilobase per Million). The False 
Discovery Rate (FDR) threshold was set to 0.1 and gene 
dispersion values were calculated by fitting a curve accord- 
ing to the DESeq "fit-only" mode. 

GO enrichment analyses 

GO enrichment analyses were conducted with the goseq 
bioconductor package version 1.10.0. As tomato data- 
bases are not yet covered by goseq built-in databases, 
transcripts lengths were retrieved with BiomaRt queries 
{Solarium lycopersicum SL2.40). An FDR cutoff of 0.1 
was used for goseq GO enrichments. Gene ontology 
terms for CC, MF and BP were similarly retrieved with 
BiomaRT queries. 

Other blolnformatic techniques 

Pie charts for the Distribution of GOSLIM BP terms were 
generated via BLAST2GO software using SL2.40 annota- 
tions. MapMan figures were obtained by first running 
the mercator tool (http://mapman.gabipd.org/web/guest/ 
mercator) with default parameters to assign MapMan bins 
to tomato transcripts. Log2 fold changes as obtained from 
DESeq output were used as MapMan input to represent 
expression changes. For clustering and heatmaps of ex- 
pression values, DESeq-normalized data were prior trans- 
formed with the Variance-stabilizing transformation (VST) 
as described by Anders and Huber [105]. Unless otherwise 
stated, further graphical outputs were generated with 
custom R and Python scripts. Custom annotations/ 
descriptions were used as specified in the text and fig- 
ures and were obtained with BLAST2GO on tomato 
SL2.40 cDNAs run with the following annotation pa- 
rameters: E-value hit filter l.E-10, Annotation cutoff 55, 
Go weight 5, Hsp-Hit coverage cutoff 20. Further com- 
plementary descriptions (ITAG2.3_descriptions) were 
obtained from ftp://ftp.solgenomics.net/tomato_genome/ 
annotation/ITAG2.3_release/. 
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Fruit mineral content measurements 

About 5 g of frozen tomato fruit pericarp from MYC 
and CONT plants were dried at 70°C for 48 hours, and 
digested in 2 ml 6 M HNO3 at 90°C for one hour. The 
digestion product was diluted to 6 ml with milliQ water 
and filtered. The P, S and K contents were determined 
in the final solution using Inductively Coupled Plasma 
Atomic Emission Spectrometry (ICP-AES) performed 
using a Liberty 100 Varian apparatus equipped with a V- 
Groove nebulizer and a Czerny-Turner monochromator 
(Department of Mineralogical and Petrological Science, 
University of Turin). The same treatments were applied 
to the control solution, which contained no sample. 
Quality control was based on the use of internal control 
samples and certified samples (Astasol-Mix from Analy- 
tica Ltd., Prague, Czech Rep.). C and N content was 
measured on fruit powder with a CHNS-O analyzer 
(Thermo Electron Corporation). 

Biochemical analyses for nutraceutical values 
Fruit extraction 

About 6 g of unfrozen fruit material were extracted on a 
magnetic stirrer for 15 min in 30 ml volume of a metha- 
nol, distilled water and formic acid (80:20:0.5) mixture. 
The extracts were double-filtered, first on blotting paper 
with a Biichner funnel connected to a vacuum pump 
and then with a 0.45 nm cellulose membrane. The total 
antioxidant activity and the phenolic content of the 
methanolic fruit extracts were immediately assayed. All 
spectrophotometric measurements were performed in 
glass cuvettes using a T80+ UV/VIS spectrometer (PG 
Instruments Ltd.). 

Total antioxidant activity 

The total antioxidant activity of the tomato fruit extracts 
was measured by testing the radical scavenging capacity 
(DPPH test) and the reducing power (FRAP assay), ac- 
cording to the methods described by Tenore et al. [106]. 

Total Phenolic content determination 

The total phenolic content was analyzed spectrophoto- 
metrically using the modified Folin-Ciocalteu method 
[107,108]. 200 \x\ of the fruit extract from MYC and 
CONT plants was mixed with 100 |il of Folin-Ciocalteu's 
phenol reagent (Sigma) in a glass test tube. Successively, 
5 ml of 7% sodium carbonate in distilled water was 
added to the tube and the mixture was kept in the dark 
at room temperature for 2 h. Absorbance was then mea- 
sured at 760 nm and results were expressed as mg Gallic 
acid equivalents (GAE) per 100 g of tomato fresh weight. 
Gallic acid (Sigma) was used to calculate the standard 
curve at 0.025, 0.05, 0.1 and 0.15 mg/100 ^il. 



Lycopene and fi-Carotene 

Quantification of Lycopene and (S-Carotene was per- 
formed spectrophotometrically according to the method 
described by Nagata and Yamashita [109]. One g of un- 
frozen fruit tissues from MYC and CONT plants was 
homogenized and extracted in 16 ml of an acetone/ 
n-hexane 40:60 solution for 15 min on a magnetic stir- 
rer. Extracts were then filtered using a 0.45 nm cellulose 
membrane and absorbance at 663, 645, 505 and 453 nm 
was measured. Contents of |3-Carotene and Lycopene 
were calculated according to the following equations: 
Lycopene (mg/100 ml) = - 0.0458 A663 + 0.3jou72 A505 - 
0.0806 A453; |3-Carotene (mg/100 ml) = 0.216 A663 - 
0.304 A505 + 0.452 A453. Results were finally expressed 
as mg/100 g fm. 

Statistical analysis 

For all the biochemical and nutraceutical measurements, 
five fruits were harvested for each MYC and CONT con- 
dition and originating from different plants. Values are 
expressed as mean ± standard deviation. Data were ana- 
lyzed with a one-way ANOVA with Tukey post-hoc or 
the non parametric test Kruskal-Wallis accepting signifi- 
cant differences at p < 0.05. Significantly different values 
are highlighted in bold. All statistical analysis were per- 
formed with SYSTAT 10 statistical software. 

Availability of supporting data 

Supporting sequence data are available in the ArrayEx- 
press database (http://www.ebi.ac.uk/arrayexpress) under 
Accession number: E-MTAB-2276. 

Additional files 



Additional file 1: Table SI. Mineral content (C, N, K, P, S) of fully ripe 
CONTROL and MYC tomato fruit Values are means ± SD (n = 5). DW, 
dry weight Data highlighted in bold are significantly different (Non 
parametric test Kruskal-Wallis p < 0.05). Table S2. Phenology and fruit 
yield of CONT and MYC plants: length of the vegetative period (expressed 
in days after transplanting, DAT), number of days from flowering to the 
first fruit, fruiting period length, average fruit number per plant, total fruit 
production, average fruit weight and total fruit yield. Data are expressed 
as mean ± standard deviation, n = 10. Statistically significant data are 
highlighted in bold according to the non-parametric Kruskal-Wallis test 
(p < 0.05). Table S4. Biochemical and nutraceutical parameters measured 
in fully ripe CONTROL and MYC tomato fruits. Data highlighted in bold 
are significantly different (one-way ANOVA with Tukey's posthoc test P < 0.05). 
Table 55. Genes and corresponding primers used in real-time RT-PCR. 

Additional file 2: Figure SI. Mycorrhization parameters of tomato 
roots inoculated with the AM fungus Funneliformis mosseae. F%: 
colonization frequency; M%: mycorrhizal colonization intensity; A%: 
arbuscule abundance; a%: arbuscule abundance in the mycorrhized 
portion. Mycorrhization parameters were determined according to 
Trouvelot method (1986). Values are means ± SD (n=5). For each plant, 
100 cm of root were measured. Figure 52. Classification of tomato fruit 
expressed genes. All expressed genes (cutoff RPKM>0.1) are represented. 
The number of transcripts assigned to tomato biological process 
(GOSLIM plant mapping) is shown in parentheses. Figure 53. 
Transcriptional changes induced by mycorrhization in tomato fruits. 
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Normalized expression mean values are plotted versus log2 fold 
changes and called DEGs (FDR<0.1) are plotted in red. Blue horizonta 
lines indicate a two-fold change in expression values for up- and down- 
regulated genes. Figure S4. Sesquiterpenoid phytoalexin biosynthesis 
pathway. Figure S5. q-PCR data for five selected genes involved in 
carotenoid metabolism on fruits from mycorrhizal and control plants. 
The y axis represents ACT value (Ct gene-Ct ubiquitin). Data are mean + 
SD for n=3. Identities of genes are: Lox (lipoxygenase TomloxC); Hdr 
(hydroxymethylbutenyl diphoshate reductase): Dxs (deoxyxylulose 
5-phosphate synthase): P^y (phytoene synthase): Ccd (carotenoid cleavage 
dioxygenase). Figure S6. Statistical analyses of q-PCR data for nine 
differentially regulated genes in three ripening stages green (G), breaker 
(B) and turning (T). The y axis represents ACT value (Ct gene-Ct ubiquitin). 
Data are mean ± SD for n = 3. One-way ANOVA (Tukey's posthoc test) 
was used for data analysis and asterisks indicate statistical significance at 
p<0.05. Identities are as follows: Solyc01g008540.2.1 (Cinnamoyl CoA 
reduaase-like protein): Solyc02g078400.2 (Allantoinase); Solyc02g076920.2.1 
(Transcription factor MYC2): Solycl0g005080.2.1 (Late elongated hypocotyl 
and circadian clock associated-1 -like): Solyc09g066410.Tl (Inorganic 
phosphate transporter): Solyc06g062540.2.1 (Putative acid phosphatase): 
Solyc03g08391 0.2.1 (Acid beta-fructofuranosidase): Solyc01g080280.2.1 
(Glutamine synthetase); SolyclOg055390.1 (Nodulin family protein). 

Additional file 3: Table S3. Genes identified as differentially regulated 
in MYC fruit in comparison with CONT fruit. Genes are ranked in 
descending order of their log2 fold change value. 
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